New insight into cataract formation — enhanced stability through mutual attraction 
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Small-angle neutron scattering experiments and molecular dynamics simulations combined with 
an application of concepts from soft matter physics to complex protein mixtures provide new insight 
into the stability of eye lens protein mixtures. Exploring this colloid-protein analogy we demonstrate 
that weak attractions between unlike proteins help to maintain lens transparency in an extremely 
sensitive and non-monotonic manner. These results not only represent an important step towards 
a better understanding of protein condensation diseases such as cataract formation, but provide 
general guidelines for tuning the stability of colloid mixtures, a topic relevant for soft matter physics 
and industrial applications. 
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It has been recognized for some time that appar- 
ently unrelated diseases including cataract, sickle-cell 
and Alzheimer's represent a broad class of 'protein con- 
densation diseases' [T]. While the study of such diseases 
has traditionally focused on detailed properties of the 
molecules involved, considerable progress has also been 
made using statistical and colloid physics. This approach 
is based on recognizing that a common feature of protein 
condensation diseases is attractive interaction between 
specific biological molecules leading to dense phases [T], 
which compromise cell and organ function. A subtle in- 
terplay between protein attractions, repulsions and en- 
tropy can lead to these condensed phases, and the differ- 
ence between health and disease can hinge on intermolec- 
ular interaction changes as small as thermal energy, kgT. 
The similar sensitivity of colloidal phase transitions sug- 
gests that colloid science tools can help understand the 
molecular origins of protein condensation diseases, and 
contribute to developing effective preventative measures. 
In a complementary fashion, globular proteins have also 
drawn great attention due to their suitability as model 
colloids [2J |21 H] ■ Proteins have proven useful for inves- 
tigating the combined effects of short-ranged attractions 
and hard and/or soft repulsions, on phase transitions, 
aggregation and cluster formation in a wide range of col- 
loidal suspensions [3J 0] . 

Cataract, clouding of the eye lens due to light scatter- 
ing and a leading cause of blindness worldwide, is an 
important protein condensation disease [TJ [6j |7j [8]. 
Statistical physics and colloid science have helped ra- 
tionalize eye lens protein solution transparency, liquid 
structure and thermodynamic properties 7 8] [TO] E] ■ 
Eye lens cells contain concentrated solutions of proteins 
called crystallins. The three major classes of mammalian 
crystallins are called a, (3, and 7 [9j [10]. The lens is 
normally highly transparent and refractive, but loss of 
transparency due to protein aggregation or phase sepa- 



ration can lead to cataract. Quantifying and understand- 
ing crystallin interactions and their impact on lens trans- 
parency is therefore a first step towards possible cataract 
prevention [BJ. 

The a-crystallins are globular, polydisperse, multi- 
subunit, 800 kDa proteins with a diameter of about 18 
nm, whose interactions are well-described with a sim- 
ple hard-sphere colloid model [12l [13] . The 7-crystallins 
are monomeric, with a molecular weight of 21 kDa and 
a diameter between 3 and 4 nm for 7B-crystallin. The 
discovery of a metastable liquid-liquid phase separation 
provided evidence for short-range attractions between 7- 
crystallins, and use of the corresponding colloid model 
has led to a quantitative description of the phase behav- 
ior nTjnanaiiEnnnis]. 

Although there are known synergetic effects in crys- 
tallin mixtures associated with the chaperone activity of 
a-crystallin, which can suppress induced aggregation of 
(3- and 7-crystallin [19] [20], little is known about the 
consequences of interactions between unlike crystallins 
at high concentrations. Recently, a-crystallin has been 
found to enhance phase separation of 7B-crystallin and 
lead to partial segregation of proteins by type in the sep- 
arated phases [21] . Therefore we have started a system- 
atic study of a- and 7B-crystallin mixtures up to concen- 
trations found within the eye lens, that combines small- 
angle neutron scattering (SANS) and molecular dynam- 
ics (MD) simulations of a coarse-grained colloidal model, 
appropriate given the hard cores of these folded, globular 
proteins. This approach allows for complementary exper- 
imental and computational investigation at the required 
large length scales and high number of proteins (~30,000 
- 60,000). 

We prepared different mixing ratios from stock solutions 
of 230 mg/mL cv-crystallin (denoted as C Q =100) and 260 
mg/mL 7B-crystallin (denoted as C Q =0) in 0.1 M sodium 
phosphate buffer in D2O at a pH of 7.1, with 20mM 
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FIG. 1: (Color online) Scattered intensities I(q) of a 230 
mg/ml (C Q = 100) and a 260 mg/mL (C a =0) solution and 
their mixtures. A, Experimental I(q) of C Q = 100 (open sym- 
bols) together with MD computer simulations (full curve). B, 
Experimental I(q) of C a ~0 (open symbols) together with MD 
computer simulations (full curve). C, /(g) of mixtures con- 
taining different amounts of the crystallin solutions shown in 
A and B. Shown are C Q = 100 (triangles), C a ~50 (diamonds), 
C a =25 (stars), C Q = 12.5 (squares), and C a ~0 (circles). Num- 
bers indicate g-values for composition dependence in Fig. [4] 



dithiothrcitol to inhibit protein oxidation and oligomer- 
ization [21]. With these precautions 7B-crystallin re- 
mained monomeric under all conditions used. We per- 
formed SANS using 1 and 2 mm Hellma quartz cells and 
varied wavelengths, sample-to-detector distances and col- 
limation lengths to cover a q-range of 0.02-3 nm" 1 . All 
experiments were performed at a temperature of 25° C, 
i.e. 10°C above the critical temperature T c of 7B- 
crystallin in the present buffer. The choice of an experi- 
mental temperature close to T c of jB-crystallin amplifies 
thermodynamic stability variations due to the presence 
of alpha crystallin, and thereby enhances sensitivity for 
discerning effective protein interaction potentials. 
We established a model for the protein interaction po- 
tentials by comparing MD computer simulations with 
the experimental scattering intensity I(q) as a function 
of the scattering vector q. Event driven MD simula- 
tions 22] were performed in cubic boxes with periodic 
boundary conditions for ^=32,000 particles for the pure 
a-crystallin solution and the mixtures, and iV=64,000 
particles for the pure 7-crystallin solution. The scat- 
tering intensities were calculated from the partial struc- 
ture factors Sij(q) (i,j=a,7) obtained from independent 
MD runs and the experimental form factors. In order 
to account for the experimental smearing, we derived a 



general resolution function and the scattered intensities 
were convoluted before comparison with SANS data [23] . 
From the three different mixing ratios (C Q =50, C a =25, 
and C Q =12, where C a is defined as 100 x [volume of 
the C a =100 solution] /[total volume of the mixture]) in- 
vestigated, C a —50 closely resembles the natural a- and 
7-crystallin concentrations found in the lens nucleus [9J. 
Therefore the model development is demonstrated on this 
sample, and its validity is tested by comparison with the 
remaining samples. 

First we determined parameters that provide a quanti- 
tative description of I(q) for the individual components. 
For a-crystallin we used a purely repulsive hard-sphere 
model and a diameter of d a = 17.6 nm [13] . which re- 
sults in perfect agreement between experimental and sim- 
ulated I(q) (Fig. 1A). For C a =0 the strongly enhanced 
intensity at low q indicates that a hard-sphere model is 
not sufficient to describe the interactions between 7B- 
crystallins (diameter c? 7 = 3.6 nm) (Fig. IB). This is a 
direct consequence of interprotein attraction and prox- 
imity to the critical concentration C c [T3J [TSJ [T7] , which 
leads to long-wavelength concentration fluctuations that 
are a potential source of lens turbidity in cataract [5]. To 
account for this we added a square-well attractive poten- 
tial with a range of 0.25 d 7 [16] and a depth w 77 = lfceT, 
where the temperature T is set to T = 0.7875 (for 
fce = 1) for all simulated mixtures. This temperature T 
has been fixed to reproduce the 7-pure case and results 
in good agreement between simulated and experimental 
I(q) in the law-q regime (q < 0.047 nm^ 1 ) relevant for 
lens transparency (FigJlB). It is clear that square well 
potentials are rather unphysical in their shape. How- 
ever, such potentials have been widely used in studies of 
colloids with short range potentials due to the fact that 
the physics of the system is scarcely dependent on the 
shape of the potential when the range is shorter than the 
diameter of the particle [23] . Figure [T] indeed demon- 
strates that simple colloid models are capable of repro- 
ducing I(q) for individual a- and 7B-crystallin solutions. 
In a next step we investigated the mixtures. Figure 
1C demonstrates the effect of adding alpha crystallin, 
from C a —0 to C Q =100. The forward scattering first in- 
creases and reaches a maximum for C a =25. For C Q =50 
it becomes highly suppressed, similar to the situation 
for pure a-crystallin. In a first attempt to simulate the 
mixtures we assumed the interactions between a- and 
7B-crystallin to be hard-sphere like. The striking dis- 
agreement between the forward scattering of the corre- 
sponding simulation with experimental data for C Q =50 
is shown in Fig. 2. The enormous increase of the simu- 
lated intensity at low q arises most likely from additional 
depletion-induced attractions known to exist in mixtures 
of hard spheres with different sizes [25]. Such an effect 
would increase T c for liquid-liquid phase separation and 
result in segregation of the proteins into large domains of 
a-crystallin-rich and 7B-crystallin-rich regions (Fig. 2). 
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FIG. 2: (Color online) Scattered intensity 7(g) of a mixture 
of a- and 7B-crystallin (C a = 50) and MD computer sim- 
ulation assuming mutual repulsion between unlike proteins. 
Left, 7(g) of C a =5Q (open symbols) together with MD sim- 
ulation of the system. Right, Snapshot from MD simulation 
shown on the left. Assuming repulsive hard sphere interac- 
tions only between a- (grey spheres) and 7B-crystallins (black 
dots) leads to a strong segregation by protein type and the 
corresponding large density fluctuations would lead to a loss 
of transparency for visible light. Also shown is a slab of the 
box for this snapshot. 



As a consequence, light scattering would have increased 
and transparency would have been lost, in contrast to the 
experimental data and the visual appearance of the sam- 
ple. Such effects are in fact manifest in a-jJ5 mixtures 
at temperatures lower than those investigated here [21] . 
However, the low forward scattering intensity 7(0) of the 
present SANS data indicates that they are suppressed un- 
der the present conditions, closer to body temperature. 
Nature must have found ways to circumvent these long- 
wavelength fluctuations and the accompanying clouding 
of the eye lens that would have been caused by hard 
sphere interactions between a- and 7B-crystallins, and 
we can speculate that an additional attraction could ex- 
ist between these unlike proteins. 

In a second attempt we thus added an attractive inter- 
action between a- and 7f3-crystallins, assuming that the 
attractive part of the potential has the same range as 
that used for the attraction between 7f3-crystallins. The 
depth of this interspecies attractive potential, u ai , was 
chosen to reproduce 7(g) and is roughly one half of the 
7B-7B-attraction. The simulated 7(g) (Fig. 3) indeed 
perfectly reproduces the SANS data throughout the en- 
tire g-range. Use of the novel attractive term leads to 
good agreement between simulations and SANS for all 
mixing ratios (Fig. 4). A snapshot of the C a =50 sim- 
ulation (Fig. 3) suggests a qualitative explanation: at- 
tractions between unlike proteins counterbalance and effi- 
ciently suppress segregation of the two proteins into large 
domains |26j . Thus transparency of concentrated crys- 



FIG. 3: (Color online) 7(g) of an a- and 7B-crystallin mixture 
{C a = 50) and simulations incorporating mutual attraction. 
Left, SANS data from Fig. 2 (open symbols) and simulations 
without (dashed curve) and with added short-range attrac- 
tions (full curve; u al — 0.55fcsT) between unlike proteins. 
Inset: 7(g* = 0.0467 nm -1 ) from simulations as a function of 
u al . The mixtures become unstable for u al < 0.5fcsT and 
> IUbT. Right, Snapshot and slab representation of the box 
for this situation as in Fig. 2. The solution is homogeneous 
and does not give rise to increased light scattering. 



tallin mixtures is maintained by introducing weak, short- 
range attractions between a- and 7-crystallins, which 
considerably decrease T c and the corresponding critical 
fluctuations at a given temperature. 
We also performed simulations in which the attraction 
between a and 7B was further increased, to above 100% 
of the 7B-7B-attraction. Such stronger attractions again 
resulted in enhanced instability (inset of Fig. 3). Thus 
the stability of these high concentration crystallin mix- 
tures depends on a-7B-attraction in a manner that is 
both extremely sensitive and non-monotonic. Such non- 
monotonic effects are a common feature of ternary liquid 
mixture phase separation [2TJ [STJ . 

There are studies suggesting that 7-crystallins inhibit 
the age-related aggregation of a-crystallins [25] . On the 
other hand, numerous investigations have suggested that 
inhibition of age- or heat-induced aggregation of 7- and 
/3-crystallins is associated with the chaperone activity of 
a-crystallins [HE HPJ US]. In contrast to these stabiliz- 
ing effects of a-7 interactions, it has recently been sug- 
gested that the observed increase in the non-covalent as- 
sociation of 7-crystallins to a-crystallins in aging bovine 
lenses might adversely affect the optical properties of the 
aged and/or cataractogenic lens [3D], since strong at- 
tractions could lead to larger aggregates and increased 
scattering of light. While such a mechanism of in- 
creased protein-protein association is a possible route 
towards cataract [3D], our experiments and simulations 
provide evidence that a weak and short-range attraction 
between a- and 7B-crystallins can potentially enhance 
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FIG. 4: (Color online) Comparison of /(g) obtained from sim- 
ulations and SANS as a function of a-crystallin mixing ra- 
tio (Cot) in the Q-7B mixtures. The predicted \ow-q portion 
of I(q) is very sensitive to alpha-gamma mutual attraction, 
unlike the higher-g portions. I(q) at different C a are from 
SANS (circles) and simulations without (diamonds) and with 
(squares) mutual attractions. Shown are data at three scat- 
tering vectors, q = 0.07, 0.2, and 0.4 nm" (the numbered 
labels refer to the arrows in Fig. [lp ). Unless explicitly plot- 
ted, error bars are equal or smaller than the symbols. 



transparency in healthy lenses. 

It is intriguing that a variety of specific molecular mech- 
anisms could potentially underly this 'stabilization due 
to mutual attraction', since a simple colloidal CV-7B at- 
traction is sufficient to rationalize our observations. Due 
in part to its small magnitude, the molecular origin of 
the inferred, net short-range attraction between a- and 
7B-crystallins remains a challenging and open question, 
as it is for short-range attractions of many globular pro- 
teins including 7-crystallin and lysozyme [21 13T]. 
It has long been recognized that uniform packing of crys- 
tallins at high concentrations, termed short-range order, 
suppresses index of refraction fluctuations and reduces 
light scattering [7.. The chemically specific origins of 
short-range order are key to understanding lens trans- 
parency. Not only do crystallins vary widely in short- 
range order properties, sensitivity to inter-crystallin in- 
teractions makes short-range order not a simple combina- 
tion of individual crystallin properties, as shown here for 
a- and 7f3-crystallins. Thus to understand transparency 
of crystallin mixtures, high concentration short-range or- 
der study, using tools of liquid-state and colloid physics, 
remains a needed complement to low concentration in- 
vestigations. 

This investigation has provided new insight into the sta- 
bility and optical properties of lens protein mixtures that 
may help understand cataract formation. The present 
results also suggest mechanisms to tune the stability of 



classical colloidal mixtures, a topic of considerable indus- 
trial importance. This demonstrates the importance of a 
colloid-oriented approach to proteins as a means of ob- 
taining insight into questions of biological, medical, as 
well as fundamental soft matter physics relevance. 
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